

NASA TECHNICAL 
MEMORANDUM 

NASA TM X-64582 


ABSOLUTE STABILITY ANALYSIS OF ATTITUDE CONTROL 
SYSTEMS FOR LARGE BOOSTERS 


By S. M. Seltzer 
Astrionics Laboratory 

D. D. Siljak 

University of Santa Clara 
Santa Clara, California 


April 2, 1971 


NASA 



George C. Marshall Space Flight Center 
Marshall Space Flight Center, Alabama 


MSFC « Fo?m 3399 {September 1968) 



TECHNICAL REPORT STANDARD TITLE PAGE 


1. REPORT NO. 2. GOVERNMENT ACCESSION NO. 

NASA TM X-64582 

3. RECIPIENT'S CATALOG NO. 

4. TITLE AND SUBTITLE 

Absolute Stability Analysis of Attitude Control 
Systems for Large Boosters 

5. REPORT DATE 

April 2, 1971 

6. PERFORMING ORGANIZATION CODE 

7. AUTHOR (S) 

S. M. Seltzer and D. D. Siljak 

8. PERFORMING ORGANIZATION REPORT # 

9. PERFORMING ORGANIZATION NAME AND ADDRESS 

George C. Marshall Space Flight Center 
Marshall Space Flight Center, Alabama 35812 

10. WORK UNIT NO. 

11. CONTRACT OR GRANT NO. 

13. TYPE OF REPORT a PERIOD COVERED 

Technical Memorandum 

12. SPONSORING AGENCY NAME AND ADDRESS 

National Aeronautics and Space Administration 
Washington, D. C. 20546 

14. SPONSORING AGENCY CODE 

15. SUPPLEMENTARY NOTES 

Prepared by Astrionics Laboratory, Science and Engineering 

16. ABSTRACT 

A method for performing an absolute stability analysis of attitude control systems for large launch 
vehicles is presented. Absolute stability of these systems is shown to be of a limited extent. The regions of 


absolute stability are computed by using a quadratic Liapunov function. The function is chosen to provide 
additional information about the exponential property of absolute stability. A system model is used to 
illustrate the method. 


17. KEY WORDS 

Booster 

Attitude control 
Stability analysis 
Nonlinear stability theory 

1.8. DISTRIBUTION STA1 

Unclassified - 

Ml dk 

-EMENT 

unlimited 

\f%/ 

19. SECURITY CLASSIF. (of this report) 

Unclassified 

20. SECURITY CLASSIF. (of this page)' 

Unclassified 

21. NO. OF PAGES 
22 

22. PRICE 

$3.00 


MSFC - Form 3292 (May 1969) 






TABLE OF CONTENTS 


Page 

SUMMARY 1 

INTRODUCTION 1 

ABSOLUTE STABILITY REGIONS 2 

LARGE BOOSTER ATTITUDE CONTROL SYSTEM DESCRIPTION 8 

APPLICATION 11 

CONCLUSION 15 

REFERENCES 16 


iii 



LIST OF ILLUSTRATIONS 


Figure Title Page 

1 . Saturn attitude control system 9 

2. Saturn navigation flow diagram ^ 

3. Simplified attitude control system 12 

4. Typical attitude control system 13 

5. Popov curve associated with equation (13) 14 

6. Nonlinearity sector 14 



TECHNICAL MEMORANDUM X- 64582 


ABSOLUTE STABILITY ANALYSIS OF ATTITUDE CONTROL 
SYSTEMS FOR LARGE BOOSTERS 

SUMMARY 


A method for performing an absolute stability analysis of attitude control systems 
for large launch vehicles is presented. Absolute stability of these systems is shown to be 
of a limited extent. The regions of absolute stability are computed by using a quadratic 
Liapunov function. The function is chosen to provide additional information about the 
exponential property of absolute stability. A system model is used to illustrate the 
method. 


INTRODUCTION 


The attitude control systems for large boosters can be analyzed in the framework 
of absolute stability. This method opens new avenues for applications of numerous 
strong results of absolute stability theory [ 1 ] to the control of large launch vehicles. 

Since the linear part of the system is not asymptotically stable, the saturation- 
type characteristic of the hydraulic actuator used for rotating a gimbaled engine violates 
the sector condition, and the absolute stability of the attitude control system is shown 
to be of a limited extent. Regions of absolute stability are estimated by using quadratic 
Liapunov functions. The Popov frequency criterion and Yakubovich matrix inequalities 
are used in selecting an appropriate Liapunov function, which provides additional infor- 
mation about the exponential properties of absolute stability. Such properties guarantee 
a desired degree of rapidity of the transient process that takes place in the system if 
perturbations occur. 

Significant advantages of the proposed method are: It is independent of the 
order of the system, and the algebraic operations involved in the computations are 
relatively simple and convenient for machine implementation. Moreover, the obtained 
results are valid not only for a particular nonlinearity but also for an entire class of 
nonlinear characteristics that satisfy certain general conditions. 

Liapunov instability criteria and frequency domain techniques Were used [2] for 
choosing the feedback control parameters in an attitude control system, which may 
result in stable operation. Conservative results were obtained so that insufficient 



information was provided for the estimation of parameters and states that assure system 
stability. In this direction, considerations of the attitude control systems for large launch 
vehicles by approximate methods [3] are more promising because they provide stronger 
“necessary conditions” in the estimation of stability regions in both the parameter and 
the state space. The approximate methods can be supplemented with the analysis of 
parametric absolute stability [ 1 ] to yield a useful addendum to the method presented 
herein. 


ABSOLUTE STABILITY REGIONS 


Attitude control systems for large boosters can be cast in the Lur’e-Postnikov 
class of nonlinear systems which are described by 1 


x(t) = Px(t) + q <t>(a), ff(t) = r^x(t) 


( 1 ) 


In equations (1), x(t) is an n vector — the state of the system; P is a constant n X n 
matrix; q and r are constant n vectors; and <t>(a) is a scalar, continuous, or discon- 
tinuous function with only isolated points of discontinuity of the first kind. 

In the absolute stability analysis of the attitude control system, the following 
aspects are essential: 

a. A necessary condition for absolute stability is that the linear part of the 
system in equations (1) described by the transfer function 


X(A) = r T (P - XI) -1 q 


( 2 ) 


be stable; i.e., the matrix P in equations (1) should be Hurwitz. However, the vehicle 
dynamics give rise to zero eigenvalues of P, and a transformation of equations ( 1 ) is 
necessary to achieve stability of the linear part. 

b. The function <p(a), which represents the nonlinearity of the system, is 
generally a saturation-type characteristic that is expected to belong to the class defined 
as 


1. Capital Roman letters denote matrices, lower case Roman letters denote vectors, 
capital Greek letters denote sets, and lower case Greek letters denote scalars. The 
letter t is used only for time, and the letter V only for a Liapunov function. Vectors 
are considered as column matrices, and superscript T denotes the transpose. The 
notation H > 0 means that H is a positive definite matrix, and I is the identity matrix. 
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<t> K : 0 < a0(a) < kct 2 , k < 00 

0(0) = 0 


( 3 ) 


Thus, the characteristic 0(a) should be situated entirely inside the sector [0, /c] of the 
a0 plane. Since the matrix P is not Hurwitz, a saturation-type characteristic violates 
the sector condition; i.e., the inequality in equation (3) and the absolute stability property 
of equations (1) are of a limited extent. Therefore, the analysis is directed towards 
estimating the finite regions of absolute stability in the state space [4] . 

c. In the application to the attitude control system, merely knowing that the 
system is stable is not sufficient. It is also important to provide information about the 
rapidity of the transient process that takes place in the system after possible perturba- 
tions. A common comparison function is the exponential function, and it will be shown 
that there exist two positive constants (p,S) independent of the initial values x(0) = x 0 , 
such that the system motion x(t) satisfies the inequality 


llx(t)ll < p llxo II e‘ 5t 


(4) 


for all x 0 , t > 0, and any 0(a) e <t> K . With this property, the system in equations ( 1 ) is 

said to be exponentially absolutely stable [ 1 ] . To achieve this kind of stability, the 
transformation mentioned in item a. must be aimed to make the matrix P + 5 1 Hurwitz. 
This, in turn, will limit the exponential property of absolute stability to a finite extent. 

We begin the absolute stability determination of the system in equations (1) by 
adding to the right-hand side the zero identity (qr^x - qa) = 0, so that equations (1) 
become 


x(t) = P tr x(t) + q0 tr (a), a(t) = r T x(t) 


(5) 


where 


P tr = p + K 5<l rT >0tr( ff ) = <H ct )- k 5 ° 


( 6 ) 


The constant «g is chosen as the least number that makes the matrix 


P 5 ~ P tr + 51 


(7) 
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a Hurwitz matrix. This is a linear problem, and any of the linear methods (e.g., parameter 
method, root locus, etc.) can be used to select the minimum value for Kg so that all the 
poles of 


XtrO) = 


X(X) 

1 +K s xO) 


( 8 ) 


are located inside the half-plane Re X < - 5 < 0. 

Once the appropriate value of Kg is found, one applies the absolute analysis to 

equations (5). The exponential property of absolute stability of equations (5) is 
guaranteed by satisfying the Popov inequality 


7t(oj) = k 1 +Rext r (" 6 + Jw) > 0 (9) 

for all real u> > 0. 

As shown by Yakubovich, the condition in inequality (9) is necessary and suffici- 
ent for the existence of a Liapunov function 

V(x) = x T Hx, H = H t (10) 

having the derivative along solutions of equations (5) as 

-V = (x^Gg x + 2x^g <f> + k -1 0 2 ) + (a - k _1 <f>) <j> + 26 V , (11) 

where 

-Cg = HP§ + Pg T H , - g = Hq + y r (12) 

and Pg is given in equation (7). 
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The matrix H > 0 satisfies the matrix inequalities 


G§ - Kgg T >0, k =£ °° 

G 5 > 0,g = 0, k = °o (13) 

Since V, -V - 26 V > 0, we can always find a sufficiently small number e > 0 
such that 


-V - 25 V > 2eV 


(14) 


By separating the variables and integrating with respect to time, we obtain 
V[x(t)] < V(x 0 )e ' 2 ( 5 +e ) t , t > 0 (15) 


Since 


llxll 2 < V(x)<m M "xll 2 06) 

where M m and are the minimum and maximum eigenvalues of H, from equation 
(15) follows the exponential property in inequality (4) of motion x(t), where 

P = (h M /h m ) 1/2 - 

To make the class <t> K of nonlinear functions 0 tr (ci) as large as possible, one 

applies the Popov graphical construction to inequality (9) and finds the maximum value 
of k for which that inequality is satisfied. 

Let us now consider a saturation-type characteristic 


0 (a) = 


fa, la I < rj 

f t? sign a, la 1 > 77 


(17) 
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where £,i? are positive constants. Because of the transformation in equations (6), the 
characteristic 0 tr O) violates the condition 

o<t>i r (o) > 0 for lal>a,a = £77 k§ -1 (1 + Kg 2 ) 1 ' 2 

which is necessary for absolute stability of equations (5). Therefore, we enlarge the class 
of nonlinear functions and require that 4>^ r (a) belongs to 

0 < a0^ r (a) < ko 2 , la I < a 

V°) = 0 a 8) 

Since Pg is Hurwitz and 0j r ( a ) e $ K a> or igi n x = 0 of the state space is 
locally exponentially stable. Therefore, a region fi§ of exponential absolute stability is 

defined as the set of all points x 0 for which the motions x(t) of equations (1) starting 
at x 0 are exponentially absolutely stable; i.e., exponentially stable for any <p lT (o) e <1> K a . 

Once the maximum k is determined such that the Popov inequality (9) is satisfied, 
and 0 tr (a) e a for some a, one is interested in finding the largest region 

S2p(J2§ 1 c £2g) defined by 

fig 1 = jx : la I < a ; V(x) < 0 j , (19) 

where the number 0 is determined from 


0 = min V(x) 
la(x)l = a 


( 20 ) 


Then, from equations (10) and (20), we obtain the conditions for 0 as 


2Hx + i>r = 0 

r T x + a = 0 (21) 
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The solution of linear equations in equation (21) for x and the constant v yields 


(3 = a 2 (r T H -1 r)' 1 (22) 


The region of exponential stability is 

fig 1 = jx : x^Hx < (3 | (23) 

and the problem is reduced to finding the matrix H > 0. 

As suggested in Reference 4, one can use the Kalman construction to determine 
a specific H from 

HPg + Pg T H = - uu T (24) 

if 7 t(co ) in equality (9) is factored as 


tt(co) = 


fl(jco) fl(- jcu) 

IPg -jwll IPg +jcoII 


and the vector u is chosen such that 


(25) 


u T (Pg -XI)- 1 q = ■ - n~ in (26) 

It is important to note that a different matrix H produces a different region fig 1 with 
respect to extent and orientation, and it is desirable to select one that provides in some 
sense the best estimate fig 1 of the region of exponential absolute stability. Then, a set 

of H matrices can be generated directly from the matrix inequalities (13), which can be 
reduced to Sylvester inequalities involving the elements of H. The best estimate 

/V/ 

fig 1 may be determined, for example, by minimizing the volume of fig 1 over the set 
of generated H matrices, which is a simple matter for quadratic regions. 
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LARGE BOOSTER ATTITUDE CONTROL SYSTEM DESCRIPTION 


A typical attitude control system for a large booster is the Saturn launch vehicle 
(Fig. 1) [5] . The major components of such a system consist of a stabilized platform, 
a launch vehicle digital computer (LVDC), an analog control computer, rate gyroscopes, 
engine actuators, and, in some instances, body-fixed control accelerometers. The 
stabilized platform is a three-gimbal device whose inner gimbal is stabilized by three 
mutually orthogonal gas-bearing gyroscopes. Mounted on the inner gimbal are three 
pendulous integrating gyroscopic accelerometers. The platform performs two functions: 
its gimbal angles provide vehicle spatial attitude, and its accelerometers detect inertial 
acceleration. The LVDC serves as the central point for all data flow and data processing; 
it solves the vehicle navigation and guidance equations, controls the sequence of flight 
events, and computes attitude error signals. As shown in detail in Figure 2, the attitude 
command signals are computed by comparing the actual attitude of the vehicle (as 
indicated by the gimbal angles of the stabilized platform) and the attitude angles required 
to guide the vehicle to its destination. The required attitude angles are calculated in the 
LVDC from the state of the vehicle as determined from stored information, stable plat- 
form accelerometer outputs, and onboard calculations indicated in Figure 2. The attitude 
error signals are converted into analog form and provided to the analog control computer 
that combines these signals with rate gyroscope signals and, in some cases (depending on 
the control law postulated), lateral accelerometer signals. From this information, it 
generates commands to the propulsion engine actuators and, when included, to the 
auxiliary propulsion system; these commands also point the thrust vector in the desired 
direction. For simplicity, it is assumed that the vehicle under consideration is rigid and 
flying above the sensible atmosphere. If rotational motion is considered in only the pitch 
plane, rotational rigid-body dynamics may be described by 


6 + c|3l — 0 


(27) 


where 


c a ig^a> (28) 

The term £ of equation (28) represents the distance between the engine gimbal point and 
the vehicle center of mass, R represents the lateral component of engine thrust, (3^ is 

the angle in the pitch plane between the thrust vector and the longitudinal axis of the 
vehicle, and I represents the vehicle’s moment of inertia about the pitch axis. The 
symbol 9 represents the pitch attitude of the vehicle. 
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INITIAL VELOCITY 



Figure 1 . Saturn attitude control system. 
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Figure 2. Saturn navigation flow diagram. 






Numerous attitude control laws have been postulated for launch vehicles and used 
successfully. One postulate may be written as 


a = a 0 'k + ai 0 + a 2 0 


(29) 


where a represents the commanded thrust vector angular rate; 'b is the measured 
attitude error; and a 0 ,ai, and a 2 are adjustable control gains. If it is assumed that the 
engine actuator produces a thrust vector angular rate proportional to the commanded 
rate until a saturation limit ( 77 ) is reached, the actuator may be modeled as equation (17), 
where <p(a) represents (3^. 

With the use of equations (27) and (29), the attitude control system (Fig. 1) may 
be simplified so that it appears as shown in Figure 3. For additional simplicity, it is 
assumed that the output of the lateral accelerometer is angular acceleration. 


APPLICATION 


Let us consider the simplified system model for the attitude control of a larger 
booster given in Figure 3. If a typical set of numerical values is chosen for c, a 0 , , and 

a 2 , such as 2, 2, 1, and 1.25, respectively, the figure may be further simplified (Fig. 4). 
This simple model, however, has all the characteristic properties of the considered class 
of systems necessary for a meaningful illustration of the outlined method of analysis. 
Since the method does not depend on the order of the system, it can be applied directly 
to more complex situations. 

For equations (1), we take 


'0 10" 


'0 ' 


"4 

0 0 1 

; q = 

0 

; r = 

2 

.0 0 0 . 


.-1 


,2.5 . 


0(ff) 


10a, la I < 0.1 
sign a, la I > 0.1 


(30) 


(31) 
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Figure 3. Simplified attitude control system. 









The transfer function in equation (2) of the linear part of the system is 


x(X) = 


2.5A* + 2A + 4 
A 3 


( 32 ) 


By applying a linear analysis, we can find that for kq 2 = all the poles of 
the transfer function in equation (8) computed as 


XtrW = 


2.5A 2 +2A + 4 
A 3 + 4.25A 2 +3.4A + 6.8 


(33) 


are located inside the half-plane Re A < -0.2. 

According to the transformation in equation (6), the system to be analyzed is 
specified by 



■0.2 

1 

0 - 


-o - 


‘4 * 

r-+ 

>-s 

II 

0 

0.2 

1 

q = 

0 

; r = 

2 


-6.8 

-3.4 

-4.25 _ 


-1 


. 2 - 5 . 


["vehicle * • I 



Figure 4. Typical attitude control system. 
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Applying the Popov graphical test to inequality (9), the Popov diagram of Figure 5 
is obtained for Xtr(^) as equation (33). As is clear from Figure 5, one can choose k = °°. 

Consequently, the system specified by equation (34) is stable exponentially in the entire 
sector [0,°°]. Because of the transformation 0 tr (a) = 0(a) - 1.7 ct, the transformed 

saturation nonlinearity violates the sector [0, °°] at la I = a = 1 1 .6 as shown in Figure 6. 
Therefore, all functions 0 tr (a) e <!>„, j j g that lie in the shaded region of Figure 6 are 
considered. 



Figure 5. Popov curve associated Figure 6 Nonlinearity sector . 

with equation (33). 


Now to determine the region J2g of exponential absolute stability for 
0 tr (a) e j j g and Xtr(^) given in equation (33), we take k = °° and factor the 

corresponding expression w(co) of inequality (9) as shown in equation (25). This 
operation yields the vector u = [-0.22 - 0.74 -2.85] ^ and the matrix 


628.46 

242.26 

184.88 

242.26 

444.02 

974.5 1 

184.88 

974.5 1 

250.65 


(35) 


as specified by equations (24) and (26). Substituting this H in equation (21), we obtain 
0 = 0.44 as defined in equation (22). Finally, we obtain an estimate of the region of 
exponential absolute stability for the system under investigation as 


m 

^ 0.2 



( 36 ) 


where H is the matrix specified in equation (35). 
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It is of interest to note that the computed region £2 q 2 in equation (36) 
corresponds to all 0 tr (a) e <1*^ j j g, which is the consequence of absolute stability. 

If we note that the transformed system is exponentially absolutely stable in the sector 
[0, °°] , which implies that the original system has the same properties of stability for 
any nonlinearity that belongs to the sector [ 1.7, °°] , we conclude that the relay charac- 
teristic 0(a) = sign a can also be used in the system in equations (1) without affecting 

the region given in equation (36). The relay characteristic, being a discontinuous 

function of the first kind, may produce sliding motions. These kinds of absolute stability 
are shown in detail in Reference 1 , and no additional analysis is necessary after the region 

£2™ is computed. 

The value of infinity was chosen for k to simplify the calculations in the example. 
However, it is reasonable to expect that better results in terms of a larger region of stability 
would accrue if a smaller value of k were chosen, as is clear from Figure 5. In other words, 
a tradeoff probably exists between the size of the region of stability and the upper limit 
of the nonlinear sector. Of course, the upper limit cannot be less than the slope of the 
nonlinearity under consideration. In practice, a factor of safety would be included by 
insuring that k is selected so that the Popov sector always includes the nonlinearity under 
study, even under the influence of expected perturbations. 


CONCLUSION 


It has been shown how the regions of exponential absolute stability in the state 
space can be computed for attitude control systems for large boosters. A quadratic 
Liapunov function was used to estimate the regions. By considering the Lur’e-Postnikov- 
type Liapunov function “the quadratic form plus an integral of the nonlinearity ,” sub- 
stantial improvements of the stability region estimation may be expected. This, however, 
would place stronger requirements either on the linear or nonlinear part of the system. 
Moreover, the numerical part of the analysis may be affected since the function is no 
longer quadratic. 

It is of interest to note that the property of exponential stability makes it 
possible to consider perturbations on the system and determine a bound on the forcing 
function that guarantees that all the solutions remain bounded inside the computed 
region. This problem is solved in Reference 4. 

Another aspect of the absolute stability analysis that is concerned with the choice 
of the vector r is not discussed herein. This problem, however, can be solved in the 
framework of parametric absolute stability as defined and developed in References 1 and 6. 
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